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Abstract 



T^ \ A diagrammatic approach based on thermal Green function to su- 

perconducting fluctuation transport is reviewed focusing on Hall and 
Nernst effects. The treatment of weak magnetic field is carefully dis- 
cussed within the linear order perturbation. 



> 

^J 1 Introduction 

This Note is the second part of the series^] and I expect that you have already 
read the first one [I]. In this second Note I shall discuss the effect of the 
magnetic field perturbationally in the weak field limit. The finite magnetic 
field shall be discussed non-perturbationally in the third Note. The symbols 
that have appeared in [I] are used here without explanations. Since the 
introduction to the series has been given in [I], let us start at once. 



*Department of Physics, Kyushu University, Fukuoka 812-8581, Japan 
1 In this Note I shall quote the first Note as [I]. 
[I] = Narikiyo: iarXiv:1112.1513l 



2 Boltzmann Transport: Relaxation- Time Ap- 
proximation 

In the relaxation-time approximation of the Boltzmann transporto the expec- 
tation value of the charge current J e and the heat current J^ in the linear 
order of the uniform electric field E and the uniform magnetic field H are 
given asj 



4^^(-^)(HxE), (!) 



and 



j0 = - 2 ^E^ 2 «p(-^)( HxE ). p) 

where both E and H are static. These describe the currents in Hall and 
Nernst effects. Although formally exact results have been obtainedjj I shall 
only discuss the results within the relaxation-time approximation in the fol- 
lowing. 



2 The Boltzmann transport is formulated in terms of E and H so that it is gauge- 
invariant. See footnote 12 in [I]. 

3 See (44) and (45) in [I]. 

4 The exact formula for J e via the Boltzmann equation is given by Kotliar, Sengupta 
and Varma: Phys. Rev. B 53, 3573 (1996) as has been discussed in §6 of [I] and it is also 
derived from the Fermi-liquid theory in [KY]. 

The exact formula for J^ via the Boltzmann equation is given by Pikulin, Hou and 
Beenakker: Phys. Rev. B 84, 035133 (2011) and it is also derived from the Fermi-liquid 
theory in [Kon]. 

In the Fermi-liquid theory the effect of scatterings beyond the relaxation-time approxi- 
mation is taken into account as the vertex correction. 

[KY] = Kohno and Yamada: Prog. Theor. Phys. 80, 623 (1988). 

[Kon] = Kontani: Phys. Rev. B 67, 014408 (2003). 



3 Quasi-particle Transport: Relaxation- Time 
Approximation 

We consider the linear response to electric field in the presence of magnetic 
field. The effect of the magnetic field is treated perturbationally in the weak- 
field limit. 

The expectation value of the charge current J 6 is expressed as the linear 
response to electric field 

J«(k,w) = ^2a^(k,u)E u (k=0,u;), (3) 

u 

where the electric field E is uniforrd_| and the k- dependence^ comes from the 
vector potential A(x) introduced as 

A(x) = Ae ikx , (4) 

with a constant vector A. The magnetic field H in the limit of k — > is 
expressed as 

H = V x A(x) = z(k x A). (5) 

The conductivity tensor is given by the Kubo formulaLJ (linear response the- 
ory) 

<V(k, w) = ^ [% u (k, uj + 15)- $^(k, iS)] , (6) 

where 

^(k,^ A )= f"dre^(T T {^(k;r)j^k = 0)}). (7) 

Jo 

Here 3 H is the charge current in the presence of the magnetic field and its 
Fourier component is given b)o 

jf(q)=^(q)-^p(q-k)A„ (8) 



5 E u (k=0,ui) should be used in (68) and (78) of [I]. 

6 We do not have to consider the k-dependence in §8 of [I] where the magnetic field is 
absent. 

7 I follow the derivation of the Hall conductivity by [FEW] . Their derivation for the 
case of nearly free electrons are extended to the case of interacting electrons by [KY] on 
the basis of the Fermi-liquid theory. These works are formulated in fully gauge- invariant 
manner. 

[FEW] = Fukuyama, Ebisawa and Wada: Prog. Theor. Phys. 42, 494 (1969). 

8 Here the diamagnetic current density is transformed by the integral 



^|d 3 a; e-^ x ^( x )^( x )A(x), 



wherqj 



j e (q) = e J2 Yl 9 ( v p + V P+q) C -'< 



t 

o- P 



cr p 



-EE P + irK^p+q- ( 9 ) 



and 



P(q) = e ^ JZ c U c p+q^- ( 10 ) 

cr P 

In this note e is chosen to be negative e < 0. It should be noted that in 
(J7J) the k-independence of j^ (k = 0) results from the coupling to the uniform 
electric field_jand the k-dependence of jjf (k; r) represents the magnetic- field 
dependence of the observed current. 

The (imaginary) time-dependence is introduced by 

J/ f(k;r) = e^ J/ f(k)e-^ (11) 

where the coupling to the magnetic fielco 

-J e (-k)-A, (12) 

is added to K in (1) of [I]. By the coupling to the vector potential the 
momentum of the electron increases by k as seen from • Thus the coupling 
leads to the electron propagator off-diagonal in the momentum. 

In the following we focus on the case of nearly free electrons where the 
electron propagator in the absence of magnetic field is given by 

G(p,i£ n ) = — — — , (13) 

It-n Sp 



with 



p(q) = e^|d 3 xe-^ x ^( x )^( x ). 



9 It should be noted that each term in j e (q) is proportional to the summation of 
incoming and outgoing velocities v p + v p+q . 

10 The uniform electric field E is expressed by the uniform vector potential Ao as 
E = zcjAo where we have put = with (j) being the scalar potential. The current 
} H '(k = 0) couples to A . 

11 It is sufficient to consider the coupling 



/>xj e (x)-A(x), 



for the discussion of the observed current linear in A. 



with 



e n = e n + — sgn(e r 
It 



(14) 



Here r is the life-time due to quasi-elastic scatterings. 







(a) 



(b) 



(c) 



(d) 



Figure 1: Diagrams with A-linear contributions: (a) The solid line in upward 
direction represents the electron propagator — G(p + k/2, ie n + iu\) and the 
broken line in downward direction represents — G(p + k/2 4— p — k/2, ie n ). 
Throughout the series of three Notes the Zeeman splitting is neglected so that 
the propagators for up-spin electrons and for down-spin electrons are degen- 
erate. The upper black circle represents a component of j^(k) where the mo- 
mentum of the electron changes from p + k/2 to p — k/2 and the lower black 
circle represents a component of j^(0) where the momentum does not change, 
(b) The broken line in upward direction is — G(p + k/2 «— p — k/2, ie n +iu>\) 
and the solid line in downward direction is — G(p — k/2, ie n ). The upper 
black circle is j^(k) and the lower black circle is j^(0). (c) The solid line in 
upward direction is — G(p + k/2, ie n + iu\) and the solid line in downward 
direction is — G(p — k/2, ie n ). The upper black circle is j^(k) and the lower 
gray circle is — (e/m)p(— k)A u where the momentum changes from p — k/2 
to p + k/2. (d) The solid line in upward direction is — G(p,ie n + iu\) and 
the solid line in downward direction is —G(p,ie n ). The upper gray circle is 
— (e/m)p(0)A M , where the momentum does not change, and the lower gray 
circle is j^(0). The integral of this diagram in terms of p is odd under the 
variable change p v — > —p u and vanishes so that (d) does not contribute to 
the conductivity. 



The Feynman diagrams for the current-current correlation function I in 



12 The procedure described in the footnote 20 of [I] can be repeated in terms of the 
off-diagonal propagator (Ti5)) . Here we put F(t,t) = (j e (k; r)j e (k = 0)). Neglecting the 



(j7j) are drawn in Fig. [T] where we only need the first order contributions in 
A to obtain the conductivity linear in H and the second order one has been 
neglected. The product of j5(k) and j„(0) leads to four kinds of terms: (i) 
JZ(k)-Jl(0), (ii) j;(k)-p(-k)A u , (iii) p(0) V#(0), and ( iv ) p(0)A,-p(-k)A u . 
In the case of (i) two current vertices can be connected by two ways as Figs. [!]• 
(a) and (b) within the linear order of A. In the cases of (ii) and (iii) one of 
the vertices is already proportional to A so that the vertices are connected by 
the electron propagator diagonal in momentum as Figs. G3-( c ) an d (d). The 



vertex correction 

F(r,r) = 2e^^ Vp+k/2 + Vp - k/2 («t_ k/2 (r)a p+k/2 (r)at,(0)a p .(0))v p , 
p p' 

is factorized as 

F(t,t) = 2e^^[<a p+k/2 (T)a p+k/2 (0))(a p _ k/2 (T)a p+k/2 (0))v p+k/2 
p 

+(a p +k/2(T)a p _ k/2 (0))(a p _ k/2 (T)a p _ k/2 (0))v p _ k/2 
within the linear order of A (See ((HI)). Since 

(a p (r)a p ,(0)) = -G(p <- p',r), <a P (r)a p (0)) = G(p, -r), 

for r > 0, 

|P _+k/2 

TO 

p-k/2 



F(r, r) = -2e ^ — [g( P + k/2, r)G(p + k/2 <- p - k/2, 
p 

+G(p + k/2 «- p - k/2, r)G(p - k/2, -r) 

771 

Introducing the Fourier transforms 



F ( r ' r ) = ^E F (^™) £ 



and 



G (P> r ) = 7j E G (P< * £ «')e~ le "' T , G(p <- p', -r) = i ^ G(p <- p', S 

n' n 

we obtain the contributions in Figs. l-(a) and (b) 

(tw m ) = -2e~ V V - [G(p + k/2, iu: m + ze„)G(p + k/2 <- p - k/2, »e n ) P + k/2 

^^ ^^ TO L TO, 

P-k/2 



F 

TO L TO 

?i p 



+G(p + k/2 <- p - k/2, iu m + ie n )G(p - k/2, ie n )- 



case (iv) is not considered here, because it is the second order contribution 
in A. 

Here G(p + k/2 ^- p — k/2, ie n ) is the Fourier transform of the off- 
diagonal propagator in momentum variable 

G(p + k/2 <- p - k/2, r) = -<T T {c p+k/2ff (r) c p _ k/2( J>, (15) 

and evaluated bvl 13 l 



G(p + k/2 <- p - k/2, ie n ) = G(p + k/2, ie n ) • ( - ev p • A) • G(p - k/2, ze n ) 

= G(p, ie n ) • ( - ev p • A) • G(p, ie n ), (16) 



whiclo is the first order perturbation in terms of (j!2p in the limit of k — >• 0. 
In order to obtain the conductivity proportional to the magnetic field 
defined in ([5]), which is linear in both A and k, we have to extract the 
k-linear contribution from the processes shown in Figs. [I]- (a), (b), (c). A 
k-linear contribution comes from the propagators diagonal in momentum 
variable p ± k/2 represented by the solid linecj It also comes from j'^(O) but 
does not from 7^(k)] 16 l 



13 If we write down the rule of Feynman diagram faithfully: 

- G(p + k/2 «- p k/2, ie n ) = [ - G(p, «£„)] • [ - (-ev p • A)] • [ - G(p, fe n )] . 

I prefer the textbook, Lifshitz and Pitaevskii: Statistical Physics Part 2 (Pergaman Press, 
Oxford, 1980), because such a faithful description is given concisely. 

14 The same relation is obtained via the gauge transformation, p — >• p eA. For 
example, in the case of free electron propagator 



G (p,ie n ) = 
is transformed into 



i£n ~(h~ " 



G (p-eA, ie n ) = G (p,ie n ) + G (p,ie n ) ■ ( p • A ) • G (p,«e„). 

15 In the case of free electrons the k-linear contribution of the free propagator is easily 
extracted as 

G (p±k/2,ie„) = G (p, ie n ) + G (p,ie n ) ■ (± — — J • G (p, ie n ). 

16 In this footnote we use diagonal propagators which appear in the right-hand side of 
(|16p . From jf,(0) vertex we obtain the contribution as 

G(p ± k/2, ie n + iu> x ) ■ j- [{p v ± y) + (p» ± y)] • G(p ± k/2, ie n ), 

7 



The k-linear part Q^Jiux) of $^(k, iuj\) is the summation of the k- 
linear contributions, $( ), $(b), $( c ), which are extracted from the processes 
shown in Fig. H]-(a), (b), (c). Since j|(k) leads to the factor (e/m)p x , jy(0) 
to (e/2m)k yi and G(p + k/2 «— p — k/2, ze n ) to C(p, z£ n ) • [— (e/m)p • A] • 
G(p,ie n ), 



$ (o) = ~^tJ2J2^{~) Pl G (P> i£ n + iux)G(p,i£ n ) 2 h y A x , 



(17) 



where 17 ! the product of the fermion-loop factor and the spin-degeneracy fac 



torj 18 ! -2, has been included. Since j* (0) leads to the facto: 



(e/2m)k. 



y 



$ 



(6) 



-2T J^ J^ - (— ) p£G(p, ze„ + ioo x ) 2 G(p, ie n )k y A a 



Since /?(— k) leads to the factor — (e 2 /m)A y , and G(p + k/2, ie n + iu\) to 
G(p,ze n + iu)\) ■ [+(l/2m)p • k] • Cr(p, ie„ + iw^), and G(p — k/2, ie n ) to 
G(p,ie n ) ■ [-(l/2m)p • k] • Gr(p,ze„), 

$ (c) = -2T^J^ [l"(~) ^(p^en + zwA^p,^)^^/ 

n p 

l/e\3 



+ o( — ) p a ,G(p,ie n + «*;A)G f (p,ie n ) fc*A/ • (19) 



which leads to a k-linear contribution 

G(p, ie„ + iu>\) ■ ( ± ^— fej/ J • G(p, «e„). 
From j® (k) vertex we obtain 



G(p + k/2, ie„ + iw A ) • — 
2m 



''" f 



G(p-k/2,fe„), 



[7 , «m 

which reduces to a k-independent contribution 

G(p,ie n + iuj\) ■ (—Pfj) -G(p,ie n ). 
\m J 

in the limit of k — > 0. 

See the footnote for the current-current correlation function in (J7|). 

17 The integrand is proportional to (p x A x +p y A y +p z A z ) -p x but the terms proportional 
to A y and A z are odd in p x and vanish by the integration over p x . In the same manner the 
integrand in (|19p is proportional to {p x k x + p y k y + p z k z ) ■ p x but the terms proportional 
to k y and k z vanish. 

18 Throughout the series of three Notes the Zeeman splitting is neglected so that the 
spin degrees of freedom only appears as the degeneracy factor 2. 

19 This factor is already proportional to k y so that we can put k = for all the propa- 
gators in the diagrams (a) and (b), because we only need the contribution linear in k y . 



Thus we obtain 



$ (o) + $ (6) + $ (c) = --(£) (k x A y - k y A x )Tj2J2ti 

n p 

x G(p, ie n + w A )G(p, ie n ) 2 - G(p, «£„ + iu x ) 2 G(p, ie r 



(20) 



If we set H = (0, 0, H) so that # = i(k x A y - k y A x ), ([20]) leads t 



H f e\3 



*?«(-») = ~ £ r ££ 



z \m 



,P a 



n p 



G'(p,2£: n + za; A )G(p,ie n ) 2 - G(p,ze n + iw A ) 2 G(p,i£ ri ) . (21) 





Figure 2: Diagrams for a fixed gauge A = (A x ,0,0): The left diagram 
corresponds to the one in Fig. [I]- (a) and the right to Fig. [T]-(b). The broken 
line represents the coupling to the magnetic field f!12p . (left) The solid line in 
upward direction is — G(p, ie n + icu\). The downward process is the product 
[—G(p,ie n )] ■ [— (— J X A X )\ ■ [— G(p, is n )]. The upper black circle is J x and 
the lower black circle is (dJ y /dp y )k y /2. (right) The solid line in downward 
direction is — G(p, ie n ). The upward process is the product [— G(p, ie n +iu>\)] • 
[— (— J X A X )] ■ [— G(p, ie n + ioj\)]. The upper black circle is J x and the lower 
black circle is —(dJ y /dp y )k y /2. Here J x = (e/m)p x and dJ y /dp y = e/m. 



If we choose the gauge A = (Aj.,0,0), we do not have to consider the 
contribution $^ of the diamagnetic current. Such a choice make the calcu- 
lation simpler] 21 ! The relevant processes leading to H-linear contribution are 



20 Eq. (2.19) in [FEW] obtained for general dispersion reduces to ([2~Tj) for isotropic 
dispersion v p = p/m. 

21 The approach by Altshuler, Khmel'nitzkii, Larkin and Lee: Phys. Rev. B 22, 5142 
(1980) considers only paramagnetic contributions $( ) and $(&) by fixing the gauge. 



9 



summarized in Fig. [2J The expressions (|17[) and (TT8|) are rewritten as 

i f) T k 

- -$ (a) = T J2 E J ^ o G ( i£ -) ( - ■/«4 B )G ! (te„)G(ze n + iw A ), (22) 



n p 



^y 2 



- -$ (fe) = T J2 J2 J *-frT ( - f ) G(ie„)G?(ie n + zwa) ( - J,A X ) G{ie n + 



n p 



zwaJ 



(23) 



where J^ = (e/m)pu and G(ie n ) = G(p,ie n ). In this case (|2ip is rewritten 



as 



^)(**) = irEE( J .) 



i <9J y 



X 



n p 

G(ie n + iw A )G(2e„) 2 - G(ie n + iu x ) 2 G(ie n ) 



(24) 



where the fermion-loop factor and the spin-degeneracy factor are moved to 
the left-hand side. Here H = —ik y A x . 

In the following I adopt a shortcmo calculation along [FEW] for 



I e (iux) = - 



-^2x(ie n )Y(ie n + iu x ), 



P 



(25) 



where X{ie n ) = [G(p,ie n )] a and Y{ie n + iu\) = [G(p, ie n + i(jJ\)} b with 
a and b being positive integers, because the calculation along §11 of [I] is 
cumbersome. By repeating the transformation which leads to (159) in [I] 



22 Using the shortcut formula (071) the summation (156) in [I] is evaluated as 
P iu + i 8)-I*(i5)=-u;f^ i (-^)G A (e)G%), 



and (167) as 



J*( W + US) I Q (*5) = - U f°° A ( - °M) e G A (e)G R {e). 



The integral 



rfe 

2ttI 



G A (e)G R {e) = iwr, 



by the residue readily leads to (165) and (171) in [I]. 



10 



and neglecting I the contributions from C\ and C4 we obtain 



/°° de r 1 

— /(e) -X A (e)Y R (e + tu x )+X A (e-tu x )Y R (e) , (26) 
-00 27R L J 

where the first integrand is the contribution along C2 and the second along 
C 3 . The contours of the integral, Ci, C 2 , C 3 and C 4 , are defined in Fig. 5 of 
[I]. To calculate the DC conductivity we only need the w-linear contribution 



P{u + i5)-I e {i5) = -uj 



—u 



de 

27ri 
de 

2vri 



f(e) 



X-(e) d -m + d -^Y«(e) 



de 



df(e) 
de 



X A (e)Y H (e 



Oe 



(27) 



where we have employed the integration by parts. By the same approxima- 
tion leading to (164) in [I], the application of ( 127|) to (12TJ) results in 



$?1)(W +«<*)- $?!)(»<$) 



~ —IUJ- 



(l) 1 



/;/ 



E 



«*l - 






c/e 
2~7ri 



G A (p,e)G*(p,e)'-G A (p,e)^(p,e) 

(28) 



where 



G K ( Pl 6) 



e-£ p -H/2r : 



G A (p,e) 



e-£ p -i/2r' 



(29) 



which is the analytic continuation of the thermal propagator ([T3"j) . The inte- 
gral over e is performed by evaluating the residue^!! so that we finally obtain 



a xy = lim lim — [^ (k, w + iS) - $^(k, «<*)] 



%* 



2 2 



(30) 



23 Here we employ the propagator (|29f and the location of its pole is explicitly known. 
The integral on C\ with G R only vanishes by choosing a closed contour in which there is 
no pole. The integral on C4 with G A only also vanishes. On the other hand, the integrals 
on Ci and C3, which contain a pair of G R ■ G A at least, are determined by the pole of 
either G R or G A . 

24 

A \G A (p, e)G R (p, ef - G A (p, e) 2 G*(p, e)l = -2r 2 . 



11 



which coincides with the result of the Boltzmann transport ([T]) taking into 
account that H = (0, 0, H) and E = (0, E, 0). 

The expectation value of the heat current J^ is expressed as the linear 
response to electric field 



jQ(k,u) = Y,a tll/ (k,u)E u (k = 0, 



(31) 



with 



cv(k, u) = - [$JUk, w + iS) - $JL(k, iS)] . 



iuj 



(32) 



As has been discussed in the above, in order to extract the contribution 
proportional to H, we only neeco the heat current J?(k) in the absence 
of the magnetic field so that only the difference between $® (k, iuj\) and 
$?y(k, iw\) is the factocj (ie n + iu\/2)/e as in the case of §11 in [I]. 



By repeating the above shortcut calculation for 

I Q (^ A ) =Aj2 fa + lJ Y)X^n)Y(lS n + i 



IU)\) 



(33) 



we obtain 
I Q {u+i5) = 



de 
2m 



f(e) 



CO 



rR, 



r+-,XV)y"u .,:) 



u> 



X A {e-u)Y R {e) 
(34) 



The cu-linear contribution becomes 

de 



I Q (u + i5)-I Q {i5) = -u 



2m 



m 



X A (e) 



dY R (e) dX A (e) 



Oe 



Oe 



-Y R (e 



X A {e)Y R {e) 



-LU 



7„ 00 2m\ de J K } 



(35) 



by the integration by parts. 

Thus the k-linear part Q^Jioox) of $^ ' (k, iu\) is obtained as 



$g ( W + iS ) - <rf [(iS) 



(1)' 



-tu- 



rn 



£ 






de 

2~ixi 



G A (p,e)G R (p,e) 2 - G A (p,e) 2 G R (p,e) 

(36) 



25 



We only need the charge current jfj(k) in the absence of the magnetic field to obtain 
<J?^„(k, w + id) proportional to H . 

26 It should be noted that the factor is proportional to the summation of incoming and 
outgoing frequencies e n + (e n + u>\). 



12 



where we have pulled out the factor —(df(e )Jd e) ■ e from the integral over e 
as in the case of §11 in [I]. Finally we obtain 2 ! 

a xy = lim lim — [$£, (k, u + iS)- $2,(k, i8)] 

~^E(-^)#^ < 37 > 

which coincides with the result of the Boltzmann transport (J2J) . 

4 GL Transport 

The linearized GL transport theory give,'' 2 ' 



a *» ~ idn?' (38) 

and 

|e| h 

axy = 4^7' (39) 



where h = 2\e\H^. The derivation of these results^] shall be discussed by 
the non-perturbational treatment of the magnetic field in the third Note. 



27 The remark in the footnote 40 of [I] also applies to (|3"7|) . 

28 See the footnote 46 in [I]. Since r 2 < in 3D as discussed in the footnote 61 of [I], 
(j3"8j) for Cooper pairs has the same sign as (|3"uT) for electrons. It is natural, because a xy 
is related to the cyclotron motion of the charged object and both electrons and Cooper 
pairs carry negative charge. By the same reason (|37p for free electrons in 3D and (|39| for 
Cooper pairs have the same sign. The same discussion also applies to a xx so that (171) 
in [I] for free electrons in 3D and (199) in [I] for Cooper pairs have the same sign. 

29 The conductivity tensor (|3"5)) and the thermo-electric tensor (|3"9")) are given in (3.52) 
and (4.36) of [3] in [I]. The contribution of the magnetization current modifies (4.36) into 
(4.38) that is identical to the result in TABLE I of [5] in [I]. 
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5 Cooper- Pair Transport: AL Process 

The calculations for electrons in §3 are translated into those for Cooper 
pairqfj straightforwardly] 31 ! The contribution of the left diagram in Fig. [3] is 

$ W = T 5Z5Z A xy^y L (^-)(- KA x )L(iu m )L(ioj m + iu x ), (40) 

m q " 

which corresponds to (j22j) and that of the right diagram is 

$ (6) =TJ2J2 ^*~df ( ~ if ) L ( iuj ^ L ( iuJm + lUJ ^ ( - KA x )L(iu m + iw A ), 

(41) 



m q a 



which corresponds to ([2"B"1) where A^ = 4eiV(0)^o9/i and L(icu m ) = L(q, zu; m ). 
Thus (1241) is translated into 

m q w 

x L(zc; m + iu x )L(iuj m ) 2 - L(iu m + iu x ) 2 L(iu m ) . (42) 

This result is rewritten as 

*fo(iw A ) = 32^(ee 2 ) 3 ^g^ e (^ A ), (43) 

q 

with 

I e (iuj\) =T^2 [L( iuJ m + iu\) 2 L(iUm) ~ L(iu m + iuj x )L(iuj m ) 2 , (44) 



30 As has been discussed in the footnote 7 of [I] , L can be identified with the propagator 
Da near T c . The order-parameter field Vt, (177) in [I], is related to the gap function A as 
* = y/zA where z = 7C(3)n/87r 2 T c 2 in 3D. 



3n 7C(3)4 e* 7C(3)n _ _e^ 

'mi& 48tt 2 T 2 9a: ~ ^J^f} qx ~ Z m* 



A% = 4eiV(0)£ 2 fe = 2e _^ 2 , l 2 ^ g a = — ^7^9, = ^— &, 



in 3D. Namely, the current vertex for \& is JJ = (e* /m*)q x in accordance with (185) in 
[I]. Here e* = 2e and m* = 2m. For example, the value of z is given in (53.23) of [FW] in 
the footnote 17 of [I]. 

31 The perturbational calculation in terms of electron propagators shall be given in the 
Supplement noticed in the footnote 63 of [I] . 
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Figure 3: AL process for a fixed gauge A = (A x ,0,0): These diagrams 
correspond to those in Fig. [2j The broken line represents the coupling to the 
magnetic field (fT2|) . (left) The wavy line in right-side is — L(q, iu m +iuj\). The 
left-side process is the product [— L(q, icu m )] ■ [—(—A X A X )] ■ [— L(q, ico m )]. The 
upper black circle is A x and the lower black circle is (dA y /dq y )k y /2. (right) 
The wavy line in left-side is — L(q, iu m ). The right-side process is the product 
[— L(q, iu) m +iu)\)] ■ [—(—A X A X )] ■ [— L(q, iu m -\-uox)\. The upper black circle 
is A x and the lower black circle is — (dA y /dq y )k y /2. Here A x = 4eN(0)^Qq x 

and dA e y /dq y = 4eiV(0)$. 



where L(iu m ) = —L(iuj m )/N(p). Analytical continuation of this discrete 
summation becomes 



where 



dx 



—n(x)l [R(+) 2 - A(-) 2 ]lmR - [R(+) - A(-)] IhlR 5 



R 



R(+) 



r] — itqX ' 
1 



.4 



Tj + ITqX ' 



Mr) 



l 



(45) 
(46) 

(47) 



7] — ir (x + to) rj + iTq [x — u>) 

Here r is complexjfj r = T\ 4- ir 2 , and T\ ^> \t 2 \. The cu-linear contribution 
is evaluated as 



p °° dr 



OX ox . 



ImR 



dR dA^ 
dx dx 



ImR' 



u> 



!°° — ( )(r-A)Ir— A— R— A—\ 
J_ 00 2iri ^ 'I & dx dx dx V' 

(48) 



32 The origin of T2 is discussed in the footnote 47 of [I] . The perturbational derivation 
of T2 shall be discussed in the Supplement noticed in the footnote 63 of [I] . 
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where ImR = (R - A)/2i and ImR 2 = (R + A)(R - A)/2i. Employing the 
high-temperature expansion] 33 ! n(x) — . T/x, 

I e (ou + i5)-I e (i5) = ujT-J, (49) 

with 

°°dxR-A(„dR A dA 8 A A dR> 



— CO 

'341 



' 2ni x I dx dx dx dx J 



Here we have picked uj 

R-A 



[t* + t q )RA, (51) 



IX 

in the integrand, because ImL R (x)/x is proportional to the Lorentz function 
in x and basic quantity. Using 

dR 9 cM „ >9 , . 

&= !T ° fl - & = -***■ < 52) 



we obtai 



r°° dr 
J = J ^(r* + r )RA{ - t R* + r*A 3 }. 

The integral is evaluated by the residua 36 ! so that 



(53) 



(r * + r ) 2 V A 2^ K J 



33 The statement between (210) and (211) in [I] is insufficient. The cut-off frequency 
u c of the fluctuation propagator, (35) in [I], is determined by the condition e = u} c tq. 
Since tq — tt/8T, uj c /T — Se/-K so that we can use the high-temperature expansion for the 
integrand with x < lo c in the limit of e — > 0. 

34 ImL R {x) = (R-A)/2i. 

35 We have used 

'"° dx(R^+A^)RA = 0. 

V ox ox J 

This is derived by the integration by parts noting that 

d ,,. jX a _ .r~dA BR- 



URA)^2RA(R^ + A^). 
ox V ox ox J 



:i(i 



da;i? 4 A = 
dziM 4 = 


2tt 1 1 


r * (i + z^ 4 ^ 

2tt 1 1 


^(i + ^V 
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Therefore 

«&£,(« + iS) - *fo(tf) = 16^e 3 u,Tr 2 £ ,^°^ 4 , (55) 

and this leads to 



q (e + e 2 q 2 ) 4 ' 



a^ = -16i7e 3 Tr 2 ^ (£+ ^ 2)4 . (56) 



q 

37 



In 2D the q-summation is performed as, 

V gogg l Tdx X l 1 (57) 



(e + £ 2 q 2 ) 4 87rdy (x + e) 4 487rcie 2 ' 

Finally we obtain ( 138|) by using ri = n/ST. 

Only the difference between <&^„(k, jwa) and ^ u (k,iu\) is the factor 
(zu; m + iuj\/2)/2e as has been discussed in §3 so that we readily obtain 

3 4(^o) 3 E^^ /0 (^), (58) 



L(iu m + iu\) 2 L(iu m ) - L(iu m + iux)L(iu m ) 

(59) 





*fi)(iu x ) = 


with 




I Q (iux) 


= T Y1 y Um + 2) 

m 



Analytical continuation of this discrete summation becomes 38 ! 

I q (lo + id) = r —n(x) (x + -) (R(+) 2 ■ ImR - R(+) ■ lmR 2 \ 

+ J" — n {x) (x - -) {iimR 2 • A(-) - ImR ■ A{-) 2 \ 

= I^(io + t5) + I 2 Q (uj + i5), (60) 



x f°° , l f°° , 1 11 

ax- — = / ax- — — el ax- 



37 

/»OC 

/ dx 

Jo (* + e) 4 7o **> + e) 3 "Jo ~> + e) 4 6e2- 

38 Previously published formulae, (35) in [Uss], the formula between (10.35) and (10.36) 
in [2] of [I] and (7) in [LNV], differ significantly from ours (|6T))) . 
[Uss] = Ussishkin: Phys. Rev. B 68, 024517 (2003). 
[LNV] = Levchenko, Norman and Varlamov: Phys. Rev. B 83, 020506 (2011). 
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where 



dx 



-n[x) x 



dx 



[R(+Y - A{-f] ImR - [R{+) - A(-)] Imi? 2 } , 

(61) 

[R(+) 2 + A(-) 2 ] ImR - [R(+) + A(-)] Imi? 2 } . 

(62) 

In the following the imaginary part of r is neglected: Tq = r or r 2 = 0. 
Employing the high-temperature expansion, n(x) '=. T/x, we obtain 39 ! 



2 ./_«, 7T 



-nix ) 



I?(u + iS)-i?(iS) = -u 

TV 



dx 



BR A dA 

2R— + 2A— 

ox ox 



ImR- 



dR dA 
dx dx 



rj 1 roc 

-u dx-^-(R + A)(R-A) 3 = 0, 

ft J -oo 2 



ImR 2 

(63) 



111 ./__ X 



| (r 2 + A 2 ) ImR - (r + a) ImR 2 } 



rp /»00 

— u / dx(-2T )R 2 A^ 

2^ J -oo 



-oo 

T 1 

"2V 



Therefore 



$g } (w + i£) - $g)(i(J) '=. 8He 2 iouT J^ 



6„2 



rt 2 



> + e 2 q 



2^2\3' 



(64) 



(65) 



In 2D the q-summation is performed ao 



£ 



4„2 



M 



(e + £ V) 3 ~ 8vrd 7 (a; + e 



dx- 



x 



1 1 



167rd e 



(66) 



39 If we put R = u{x) + iv(x) and A = u(x) — ii>(a:), u(x) is even: u(— x) = u(x) and 
«(&) is odd: «(— x) = — v(x) in a;. Namely R + A is even and i? — A is odd. Therefore the 
integrand (R + A)(R — A) 3 is odd so that the integral in ((63)) vanishes. 

f°° , 2 .2 Til 

./-oo 2 T 7^ 



10 



dx 



{x + ef 



1 f°,J X - 11 

dX Jx~Te)^' € J dX Jx~Tef-27- 
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Finally we obtain ( 159"j) via 



a 



■xy 



e 2 1 

T h H 4 



(67) 



The calculations in the absence of the magnetic field are also performed 
in the same manner. The integral (207) in [I] with real Tq is evaluated aa 41 l 



I e (u + i5)-r(i5) = 



T 



N(0) 2 n 
1 T 

N(0) 2 n 
T 



UJ 



dx 



dR 

dx 



dA 

dx 



R-A 



IUJTq 



2ix 
dx(R 2 + A 2 )RA 



2ivW^V' 



The integral (216) in [I] with complex r is evaluated aa_f 



I Q (uj + i5)-I Q (i5) = 
+ 
to give (220) in [I]. 



T 



UJ 



N(0) 2 vr 

1 Tu 
N(0) 2 n2 



dx 



-oo 

'OO 



dR 

.dx 



dA 
dx 



R-A 
2% 



dx(R - A) 



R-A 

2ix 



(6* 



(69) 



6 Remarks 

The sections of Exercise and Acknowledgements are common to [I] so 
that I do not repeat here. Some typographic errors in [I] are listed below. 



dxR A A 



dxRA i 



TT 1 1 
4 T 7? 3 
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Noting that 



d , A\ 2 -o(n _ A\(— - —\ 

dx V dx dx J' 



the integral in the first line of ([5!)| is shown to vanish by integration by parts. The integral 
in the second line is evaluated using ([5Tj) and 



dxR 2 A = 



dxRA 2 



2tt 


1 


1 


T 


( 1 + t) 2 


rf 2 


2tt 


1 


1 



ro^ + J)^ 2 ' 
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• The thermodynamic relation in p. 11 should be dQ = dU — fidN. 

• The electric field in (68) and (78) should be E u (k—0,uj). 

• The factor e in the first line of (73) should be e/m. 

• In Fig. 5 the subscript A for the lower cut C\ has dropped by the 
font-error at arXiv. 

• In the footnote 47 the first term in the right-hand side of the complex 
GL equation should be (1 + ic )\I/. 
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